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A light vector boson, Zd, associated with a "dark sector" U{l)d gauge group has been introduced 
to explain certain astrophysical observations as well as low energy laboratory anomalies. In such 
models, the Higgs boson may decay into X ~\-Zd, where X = Z, Zd or 7. Here, we provide estimates 
of those decay rates as functions of the Zd coupling through either mass-mixing [e.g. via an enlarged 
Higgs mechanism) or through heavy new fermion loops and examine the implied LHC phenomenol- 
ogy. Our studies focus on the higher m^^ case, > several GeV, where the rates are potentially 
measurable at the LHC, for interesting regions of parameter spaces, at a level complementary to 
low energy experimental searches for the Zd- We also show how measurement of the Zd polarization 
(longitudinal versus transverse) can be used to distinguish the physics underlying these rare decays. 



I. INTRODUCTION 

The discovery of a 125 — 126 GeV scalar state at the 
Large Hadron Collider (LHC) [l|, 0) appears to have pro- 
vided the last missing ingredient, namely the Higgs bo- 
son, of the very successful Standard Model (SM) of parti- 
cle physics. Given its important connection to generation 
of elementary particle masses, the new state, henceforth 
referred to as the Higgs and denoted by iJ, is poten- 
tially a good place to look for physics beyond the SM. 
For example, _ff — >■ 77 can deviate from the SM expec- 
tation [3, IJ] if new weak scale states contribute to the 
loop processes that mediate the decay. Early hints of 
such a deviation may already be present. However, fur- 
ther analysis and more data are necessary to address this 
question, as the current experimental situation does not 
allow one to draw definite conclusions |5,] . 

In any event, the Higgs could provide other surprises 
with more data, since there are good reasons to expect 
physics beyond the SM is required to explain Nature. 
Chief among them is the presence of dark matter (DM), 
accounting for more than 80% of the matter in the Uni- 
verse Q. 

On general grounds, one may expect that the DM is 
part of a larger particle sector whose interactions with 
the visible matter (SM) are not completely decoupled. 
Such a point of view has been adopted in explaining some 
astrophysical data that could be interpreted as DM sig- 
nals. For examples, see Refs. |7J-y] for an explanation 
of the 511 keV gamma-rays observed by Integral/SPI 
[10| and Ref. [il[ for an explanation of the electron and 
positron excesses observed by ATIC [1^ and PAMELA 
|13| . Independent of DM considerations, there are also 
other possible clues that point to the presence of new 
physics, in particular the 3.6cr discrepancy between the 
measured value of the niuon anomalous magnetic mo- 
ment g^i — 2 and its predicted value in the SM [^, [Ij, [I^ . 
This and other tentative, less significant hints of beyond 
the SM physics may also potentially originate from the 
same "dark sector" that includes DM. 

In light of the above considerations, we examine here 



how aU{l)d gauge interaction in the "dark" or "hidden" 
sector may manifest itself in Higgs decays. In particular, 
we consider the possibility that this Abelian symmetry 
is broken, giving rise to a relatively light vector boson 
Zd with mass mz^ < 10 GeV. We further assume that 
Zd, being a hidden sector state, does not couple to any 
of the SM particles including the Higgs directly, i.e. SM 
particles do not carry dark charges. We do allow for 
the possibility of particles in extensions of SM, such as 
a second Higgs doublet or new heavy leptons, to carry 
dark charges, leading to indirect couplings, via Z — Zd 
mixing or loop-induced interactions. This allows us to 
assume that the properties of H are close to those of the 
SM Higgs, as the current experimental evidence seems to 
suggest. For simplicity we will ignore the possibility of 
H mixing with other scalars unless specified otherwise. 
We will consider decays of the type 



H^xZd {with x^z,Zdn), 



(1) 



and examine the prospects for detecting such signals at 
the LHC with emphasis on the H — >■ ZZd decay mode. 
Other possibilities will only be briefly commented on. 



II. FORMALISM AND PHENOMENOLOGY 

Quite generally, there are two classes of operators that 
contribute to the H decays of interest in our work: (A) 
dimension-3 and (B) dimension-5 operators. 

A. Dimension-3 operators 



We have 



Oa,x = CA,xHXfj,Zj^ , 



(2) 



^ We will not separately list operators that involve derivatives of 
the Higgs llal . as they can be recast, using equations of motion, 
into one of the forms considered below. 
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FIG. 1: Processes contributing to interactions for (a) case (A) 
of diniension-3 operator (tree-level mixing) and (b) case (B) 
of dimension-5 operator (loop-induced decay through a new 
fermion F). 



with ca^x a coefficient that has mass dimension +1. By 
gauge invariance X ^ 7 in this case. This kind of oper- 
ator typicaUy arises from mass mixing between the SM 
Z and Z^, or by Higgs mixing with another scalar, such 
as a second Higgs doublet, that carries C/(l)d charge and 
gets a vacuum expectation value (vev) . [See Fig. [1] (aV| 
For a discussion of the latter possibility see, e.5., Ref. [I4I 
where Higgs mixing with a SM singlet scalar charged un- 
der U(\)d has been considered. Whatever its origin, for 
processes mediated by the operator Oa,x, the boosted 
Zd final states will be primarily longitudinal. 

For operators of type Oa,x in Eq. ([5]), we wiU focus on 
X = Z . Such interactions are typically associated with 
mixing. For example, the mass term for Z-Z^ mixing can 
be parametrized as e^m^ZZc;, with 



ez 



mz 



S, 



(3) 



where (5 is a model dependent parameter. (We will as- 
sume that Zd is light compared to Z in this work.) The 
vector boson Zd couples to the weak neutral current, like 
the SM Z boson, with a coupling suppressed by ez- This 
coup ling can then provide a new source of parity viola- 
tion [18| and couplings of Zd to non-conserved currents. 
At energies E > tti^^, the longitudinal Zd has enhanced 
couplings of order (E/mz^)sz which dominate its weak 
scale processes, as implied by the Goldstone boson equiv- 
alence theorem [l^. Therefore, we expect the Higgs de- 
cay H — )■ ZZd mediated by Oa,x to be dominate d by 
longitudinal vector bosons. As discussed in Refs. [l8l.[20|. 
this scenario can be constrained by various experiments, 
both at low and high energies. (See Ref. [2l| for a re- 
cent review of the relevant low energy parity violation 
experiments.) 

Roughly speaking, rare B and K decays suggest S^ < 
10~^ for mzj ^ 5 GeV (a scale set by the B meson 
mass) while the good agreement between precision Z 
pole property measurements and SM predictions imply 
S^ < few X 10"'' for all mzj [IS]- For comparison, 
the good agreement between early H decay data at the 
LHC already suggests 6'^ < 10""' for some values of mz^ 
[18] which is competitive with precision Z pole studies 



and the (somewhat model dependent) rare meson decays. 
Hence, continued future searches for H ^ XZd hold the 
promise of providing a sensitive, unique probe oi Z — Zd 
mixing, particularly in the higher mz^ region above a few 
GeV where rare flavor decays are not applicable. 



B. Dimension-5 operators 

These operators can be CP even or odd. The CP even 
interaction has the form 



Oi 



,x 



CB.xHXnpZ 



M^'^d ' 



(4) 



where X^^„ ~ dfj_X^—d^X^ is the field strength associated 
with X. The CP odd interaction can be written as 



O 



cb.x 



B,X 



Sfii^paHX^'^Zd'"^ 



(5) 



Here, cb,x and cb,x have mass dimension —1 and are 
model dependent. Generically, these interactions result 
from 1-loop processes [see Fig. [T] (b)]. For example, 
in models with vector-like fermions carrying SU{2) x 
U{1)y X U{l)d quantum numbers such as the heavy lep- 
ton model of Ref. ^^ , these operators can naturally oc- 
cur. The CP odd interaction in Eq. ([5]) arises if the new 
fermions have complex couplings to the Higgs that intro- 
duce CP violating physical phases |23{ . 

One-loop effects from the aforementioned vector 
fermions typically result in kinetic mixing between U{l)d 
and U{1)y [24|, parametrized by 



2 cos 61 



w 



^fiW^d ■ 



(6) 



where 9w is the weak mixing angle and B^,j is the hy- 
percharge field strength tensor. For a loop-generated e, 
we would generally expect e ~ e5d/(167r^) '^ O(10~'^) 
assuming unit charges; e is the electromagnetic coupling 
constant and gd is the U{\)d gauge coupling. The vec- 
tor Zd can couple to the electromagnetic current like a 
photon with a coupling suppressed by e. 

A light Zd with small coupling to the electromagnetic 
current can contribute to the 3^ — 2. The relevant dia- 
gram, explicitly showing the ^-Zd mixing through heavy 




FIG. 2: Zd contribution to the g^ — 2. The new fermion (F) 
loop-induced ^-Zd mixing is expUcitly shown. 




500 1000 




50 100 
»!zci [MeV] 
(b) 



500 1000 



FIG. 3: (a) Shaded regions of £ — mz^ parameter space ruled out by various experimental constraints. The green band shows 
the parameter region (90% CL) that can explain the 3.6cr discrepancy in g^ — 2 [23]. (b) The same parameter space showing 
the sensitivity of future direct light Zd gauge boson searches. Such search constraints generally assume Br(Zd — >■ i^£^) ~ 1. 



fermion (F) loops, is given in Fig. [21 A similar loop- 
induced vertex amplitude is illustrated in Fig. [1] (b) for 
the Higgs decay into Z^. 

The Zd contribution to g^ — 2 can explain the current 
3.6tT deviation of the muon anomalous magnetic moment 
from the SM prediction |25l . |26| in a parameter region 
indicated by the green band in the e^ — rnz^ parame- 
ter space in Fig. [3| (a). As the figure shows, roughly 
£2 _ 20-6 - 10-5 and niz^ ~ 20 - 50 MeV is left, after 
various experimental constraints with the most recently 
updated values are imp osed: electron g — 2 [22, [23| , beam 
dump experiments [2^ [23], meson decays into Zd {T de- 
cays at BaB ar t28il . ^ decays at KLOE [30|, tt" decays at 



SINDRUM |3ll32| and WASA-at-COSY [33]), and fixed 
target experiments (MAMI (H], APEX [35]). The req- 
uisite value of e for the t;^ — 2 explanation can naturally 
arise from 1-loop diagrams ii gd ~ e. Figure [3|(b) shows 
the expected sensitivities of future Zd searches beyond 
existing bounds in the same parameter space [36|. 

As we discussed, the low mass region {mz^ < 1 GeV) is 
particularly well motivated by the g^ — 2 discrepancy but 
it is being thoroughly explored by direct production at 
electron facilities (such as the one at Jefferson Lab in the 
US, and the one at Mainz in Germany), as well as rare 
meson decays. In this paper, we will concentrate on the 
higher mass region {mz^ ^ 1 GeV) which is potentially 
accessible in Higgs decays. 

Assuming that the couplings in Eqs. ^ and ([S|) are 
induced by loops of vector-like fermions, with SU{2) x 
U{1)y preserving masses mp, we expect cb,x, cb.x — > 



as nip — >■ (X). The masses of these fermions are shifted af- 
ter electroweak symmetry breaking, due to their Yukawa 
interactions with the Higgs. For mp ~ few x 100 GeV 
and Yukawa couplings yp ^ 1 the lightest new fermions 
can have masses of 0{niz) and [for 0(1) phases] we get 



\CB,X\ 



\CB.X\ 



9wgd 

IGtt^ mz 



(7) 



with g^ a typical electroweak coupling constant. Such 
vector-like fermions [but not necessarily charged under 
U{l)d\ have been recently proposed as a natural solution 
to the larger than expected diphoton rate in the early 
Higgs data at the LHC (for some examples, see Refs. (371 - 

US)- ^ 

As discussed in Ref. [22|, with the additional assump- 
tion of being charged under U{l)d, the aforementioned 
fermions are also well motivated as mediators of the 
U{1)y and U{l)d kinetic mixing. In this framework, 
we can estimate the H — >■ "^Zd and H — > ZdZd rates, 
using ^-Zd mixing. For gd ~ e, as motivated by the 
above discussion of g^ — 2, it is estimated that Bt{H — >■ 
ZdZd) < Bt{H -^ -fZd) ~ 0.1Bt{H -^ 77) if the vector- 
like fermions are to explain the Higgs to diphoton decay 
rate at about 1.5 times the SM prediction^. 



^ At leading order, W loops do not contribute to Ob,x and Ob,x 
while their contribution is dominant over the fermion loop con- 



tributions in the SM H — > 77, Zj processes. The H 
rates would then likely be comparatively suppressed. 



XZd 



Given the preceding discussion, for Br(Zd -^ £^£~) > 
Bt{Z — > £~^£~), one could expect signals of iJ — > XZ^ 
from the interactions in Eqs. (jlj and ([5]) at the LHC, 
with current or near future levels of statistics. We also 
note that the final state Zd from decays mediated by 
Ob,x and Ob,x would be primarily transversely polar- 
ized. The Zd can mimic the promptly converted photon 
(7 — >■ e^e~) although the small nonzero mass and a pro- 
duction vertex near the beam can be used to distinguish 
the new Zd events. 



III. PARAMETERIZATION AND SETUP 

Models that could give rise to operators Oa,x, Ob,x, 
and Ob X have recently been studied, for example, in 
Refs. [i^ [22, [23, [4^ . Here, we briefly explain our param- 
eterization and set up our notation. The general decay 
widths for each case can be found in our Appendix. 

Class ( A) : Regardless of the origin of Oa,x (be it mix- 
ing between scalars or vectors), following Z — Zd mixing 
we will parameterize the interaction in Eq. ([2|) by 



H— > Z Z -^ e e' |i |i' 



Ca,X 



COS 6w 



£zmz , 



(8) 



where g is the SU{2)l couphng constant. 

Class (B): In analogy with ca,x, "we parameterize the 
strength of interactions from Ob,x and Ob by 



and 



cb.x 



cb.x 



2 cos 6* 



w 



-{nx/mz) 



g_ 

2 cos 61 



w 



-{kx/mz) , 



(9) 



(10) 



respectively, where kx and kx are dimensionless model 
dependent constants. The normalizations have been cho- 
sen for ease of notation in our later results. 

As previously pointed out, the interactions mediated 
by class (A) and (B) operators result in final states that 
are dominated by the longitudinal and transverse polar- 
izations of the vector bosons, respectively. This effect 
can be used as a diagnostic probe of the underlying mi- 
croscopic process to determine whether they result from 
class (A) or class (B). The operators Ob.x and Ob,x 
can be further disentangled if the angular distribution of 
the final state particles are considered. For the Higgs to 
diphoton decay, for example, see Ref. |23| . 

The mass range we consider in our LHC study is pri- 
marily m^d ~ 5 — 10 GeV which is complementary to 
searches of lighter Zd in low energy experiments; see 
Fig. [3] Our examination of Higgs decays to light vec- 
tor bosons at the LHC is also complementary to studies 
based on heavy (TeV scale) Z' scenarios (for example, 
see Refs. |49|,[50[). Before closing this section, we would 
like to point out that the size of 'Br{Zd — > £'^£~), £ = e,fi, 
is model dependent. For example, if kinetic mixing be- 
tween U{l)d and U{1)y {s 7^ 0) dominates, one might 
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FIG. 4: Transverse momentum distribution of the hardest 
(solid) and softest (dashed) leptons identified as originating 
from the Z (black) and Zd (red/grey) for the e^e~ /i^ jj.~ final 
state. We use the parameter choices of Eqs. (|12l I13[l . No 
energy smearing or cuts have been applied. 



expect Br{Zd ^ £+£^) ~ 0.3 for the 5-10 GeV mass 
range while for Z — Zd mixing dominance it could be 
smaller [18'|. Because of that arbitrariness, we will give 
our results in terms of 6^ x Br {Zd — >■ £~^£~). 



IV. NUMERICAL ANALYSIS 

We now present the results of our collider simulation. 
We study the observability of 



PP 



H -^ ZZd -^ £t£iit^2 



(11) 



at the y/s = 14 TeV LHC, where £1 and £2 are electrons 
or muons. 

The decay H — >■ ZZd is an interesting channel to study 
as it can arise in both class (A) and class (B) type in- 
teractions. The H — >■ jZd mode is not relevant for class 
(A), since gauge invariance does not allow it. While the 
H -^ ZdZd decay may also occur for both classes of inter- 
actions, it will be doubly suppressed if it results only from 
Z-Zd mass mixing. The decay widths for all three modes 
are presented in the Appendix, using the parametrization 
of the previous section (ez, nx, kx)- 

The region of interest is a relatively light Zd, and so, 
unless otherwise noted, we use a Zd and Higgs mass of 



mz^ 



5 GeV and tuh = 125 GeV. 



fl2) 



In our analysis we find that the kinematic cuts affect 
events from operators of class (A) and (B) somewhat 
differently. For simplicity, the results of the numerical 
simulation are presented in detail with the coupling con- 
stants 



S^ X Br{Zd -^ £'^t 



IQ- 



and Kz — kz — 0. (13) 



More general cases can be obtained by rescaling the over- 
all coupling and taking into account how the kinematic 
cuts effect the two classes of events, as is detailed at the 
end of our analysis. 

Using the results of Ref. 18], the above choice corre- 
sponds to Br(iJ -^ ZZd) « 2 X 10"'*/Br(Zrf -^ £+(.-). As 
mentioned before, when kinetic mixing dominates, the 
branching ratio of Zd into leptons is typically given by 
Br(Zd -^ e+e") ~ Br(Zd -^ M^M") - 0.15. 

We now discuss how to isolate our signal from back- 
grounds. All signal and backgrounds were simulated us- 
ing MadGra^h 5 [U with the CTEQ6L parton distri- 
bution set [53 ■ The operators in Eqs. (|2I4I5P with the 
parameterization in Eqs. (l8]- lip[) w ere implemented into 
MadGraph 5 using FeynRules |53J |. 

To isolate the signal, it will be useful to implement full 
event reconstruction. Hence, we need to identify which 
opposite sign, same flavor lepton pair originates from the 
Z and which originates from the Zd- Since both the Z 
and Zd have narrow widths, we expect one pair of lep- 
tons to reconstruct the Z mass and the other pair the 
Zd mass. Since the Zd mass is not known a priori^ the 
Z boson mass is utilized to identify the decay products. 
We calculate the invariant mass of all opposite sign, same 
flavor lepton pairs and identify the pair that closest re- 
constructs the Z pole as originating from the Z boson. 
The remaining pair is identified as originating from the 
Zd- We will sec that the invariant mass distribution of 
the lepton pair identified as originating from the Zd will 
permit a measurement of mz^ ■ 

Detector resolution effects are simulated via Gaussian 
energy smearing applied to leptons and jets with a stan- 
dard deviation parameterized as: 



E 



E 



(14) 



where (j{E) is the energy resolution at energy E, © rep- 
resents addition in quadrature, and all energies are mea- 
sured in GeV. For leptons (jets) we take the ATLAS val- 
ues [El a = 10% (50%) and b = 0.7% (3%). The energy 
resolution of electrons and muons have different depen- 
dencies on EM calorimetry and charged particle tracking. 
However, we use the uniform values for EM calorimetry 
energy resolution for all final state leptons, which, for the 
energies under consideration, is more conservative than 
tracking capabilities. 

For events to be triggered on, they must pass mini- 
mum acceptance cuts on the rapidities, 77, and transverse 
momenta, pT, of the final state particles. Typically, for 
leptons a pT cut around 20 GeV is used. However, for 
our signal, two of the leptons originate from a low mass 
resonance and are expected to have low transverse mo- 
menta. In Fig. 21 we show the transverse momentum 
distributions of the hardest (solid) and softest (dashed) 
leptons identified as originating from the Z (black) and 
Zd (red/gray) for the e~^e~iJ,~^^~ final state. Distribu- 
tions for the 26+26" and 2/x"'"2y^~ final states are similar. 
The leptons from the Z decay are typically harder than 



those from the Zd- This can be understood by noting 
that the momentum of the Z and Zd in the Higgs rest 
frame is \p\ « 30 GeV. Hence, the energy of the Z is 
dominated by the Z mass, mz, while the energy of a 
light Zd is dominated by \p\. The decay products of the 
Z then have higher transverse momentum than the decay 
products of the Zd and typically peak near mz/2. Based 
on these considerations, we apply transverse momentum 
and rapidity cuts on all final state leptons: 



p^ > 4 GeV and \ri'^\ < 2.5. 



(15) 



To trigger on an event, ATLAS and CMS typically re- 
quire that at least one final state particle have pr larger 
than that required by Eq. ^. We follow the ATLAS 
triggers [l] on iJ — > ZZ* -^ 2i^2i~ events and require 
that one lepton passes a transverse momentum thresh- 
old of 24 GeV or two leptons have a minimum transverse 
momentum of 13 GeV each. 

Signal event characteristics are exploited to develop 
additional cuts: 

• To resolve the four leptons of the signal, any pair 
of leptons is required to satisfy the isolation cut 



Ai? ^ \/(A?7)2 -H (A0)2 > 0.3, 



(16) 



where A77 and A0 are the rapidity and azimuthal 
angle differences between the two leptons under 
consideration. 

• Since signal events originate from a Higgs boson, 
the invariant mass of the four leptons should closely 
reconstruct the Higgs mass. Using typical mass 
resolutions [l| for a 125 GeV Higgs, the four lepton 
invariant mass, 7714^, is required to satisfy 



|m4£ — m//| < 2 GeV. 



(17) 



• Since the Z has a relatively narrow width, the lep- 
ton pair identified as originating from the Z should 
reconstruct mz- Hence, we require 



mz\ < 15 GeV, 



(18) 



where m^£'^ is the invariant mass of the recon- 
structed Z boson. 

The main backgrounds to our signal consist of the ir- 
reducible backgrounds 



H^ZZ*, Z7*, ZZ, Z{-^Ae), 



(19) 



the reducible backgrounds Zjj, where the Z decays lep- 
tonically and the jets fake leptons, and the ti decay chain 



tt 



cr){b^c£+)e+e- +^T, 



(20) 



where ^^ is missing transverse momentum. We assume 
that a jet fakes an electron or muon 0.1% of the time [3| 
and use the branching ratio Br(6 — )■ c £D) ~ Br(i?° — )■ 
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FIG. 5: Distribution of the reconstructed Zd invariant mass, 
for the e^e~ fi'^ fi~ final state. The parameters of Eqs. H12I 
I18[> were applied. 



mTz^, for both major background (shown cumulatively) and signal 
[T3)l were used. Energy smearing, triggers, and cuts in Eqs. (fTSl- 



TABLE I: Signal and background cross sections (in fb) with consecutive cuts, and signal to back; 
for the e'^ e~ fi'^ n^ , 2/i+2/i^, and 2e'^2e~ final states at ^/s = 14 TeV. The parameters in Eqs 



Channel 


e+e- 


U'^fJ- 


2/^+2/i- 


2 


2+2e- 


a(fb) 


Signal 


Background 


Signal 


Background 


Signal 


Background 


No cuts and no energy smearing 


0.10 




0.051 




0.051 




Basic cuts (1151) + Triaeer + Isol. (Il6l) 


0.049 


67 


0.024 


26 


0.024 


26 


+ m« (mi) + mT HSI + mF; (EH 


0.043 


0.030 


0.022 


0.017 


0.022 


0.014 


S/B 


1.5 


1.3 


1.5 1 



jround ratio S/B after all cuts 
T2l[T3l) were used. 



£D + anything) = 0.10 [y]. Also, the analysis of the tt 
background only takes into consideration the four lep- 
tons in the final state, ignoring the missing transverse 
momentum and extra hadronic activity from the b de- 
cays. 

The first row of cross sections in Table H] shows the sig- 
nal rate before energy smearing, while in the second row 
the signal and background rates are given after energy 
smearing, triggers, and cuts in Eqs. (|15ll6p . Background 
cross sections before minimum cuts are not shown since 
soft and coUinear singularities in some channels do not 
allow for a reliable estimate. After these cuts, the dom- 
inant backgrounds are ti and Z , making up ~ 50% and 
^ 28% of the e^e~ fX^ ii~ background and ^ 32% and 
~ 38% of the same (lepton) flavor background, respec- 
tively. The next largest background is ZZ , contributing 
~ 12% to e"'"e~^"'"/i~ and ~ 26% to same (lepton) flavor 
backgrounds. Unlike the other backgrounds, which re- 
quire that at least one same flavor lepton pair originate 
from a Z , any flavor lepton can originate from a 6 or W^ in 
the top decay. Hence, the different percent contributions 
to e^e~ fjL^ pL~ and same flavor final states for different 
backgrounds are due to the combinatorics of the tt de- 



cay into four leptons. The tt, Z, and ZZ backgrounds 
are mostly eliminated by the mass cuts in Eqs. (|17ll8p . 
with the large majority of the leftover background com- 
ing from Z"f* andH -^ ZZ*. [See Fig. (P.] 



Zd Resonance Peak 



In Fig.[5l we show the reconstructed Zd invariant mass, 
m^£'^, distribution for major backgrounds (shown cumu- 
latively) and our signal for the e~^e~ ^'^fj,~ signal after 
energy smearing, triggers, and cuts in Eqs. (fTSJ - [T51) . 
The distribution is similar for the 26+26" and 2/i,+ 2/z^ 
signals. As can be seen, the signal clearly stands out 
above background at an invariant mass of mz^ = 5 GeV, 
and most backgrounds, except Zj* and H — >■ ZZ*, are 
negligible after the mass cuts in Eqs. (|17ll8p . The dilep- 
ton distributions from H -> ZZd and H — >■ ZZ* in Fig. [5] 
show agreement with the results of Ref. [l^ which was 
based on Ref. [55j . 

The sharp fall in the invariant mass of the background 
below 5 GeV can be understood by noting the invariant 



mass of two massless leptons can be written as 
mi2 = 2£;i^2(l - cos6'i2), 



(21) 



where E12 and O12 are the energies and angular sepa- 
ration of the two leptons, respectively. Cuts placing a 
minimum on the energy and the angular separation of 
the two leptons can be reinterpreted as a minimum on 
the dilepton invariant mass. The cuts in Eqs. (|15|16|) 
effectively place a minimum on m'^^'^, as seen in Fig. [S] 
Hence, to probe Zd masses below '-^ 4 — 5 GeV at the 
LHC, either the isolation or transverse momentum cuts 
need to be relaxed. 

Motivated by Fig. [SJ we place the cut on the Zd recon- 
structed mass: 



mz^ < O.l mz^ 



(22) 



To illustrate the efficiency of these cuts. Table U lists 
the cross sections for both signal and background bro- 
ken down by dilepton final state signals. The second row 
shows the cross sections after energy smearing, triggers, 
and cuts in Eqs. (|15ll6p . while the third row includes the 
additional invariant mass cuts of Eqs. (I17I18I22I) . In the 
last row we show the signal to background ratio, S/B, 
after all cuts. As can be clearly seen, the invariant mass 
cuts leave the signal cross section mostly intact while 
severely suppressing the backgrounds, and the S/B ratio 
is well under control. 

We also explore how sensitive the LHC is to the pa- 
rameters mza and S^ x Br{Zd -^ i^£~). In Table HIl we 
illustrate the luminosities needed at -/s = 14 TeV for 
2(7 exclusion, 3ct observation, and 5cr discovery for both 
mza = 5 GeV and mz^ = 10 GeV. Significances have 
been calculated as S/\/S + B. All other parameters are 
the same as in Eqs. ([121 [T3|) . For a slightly heavier Zd the 
LHC is found to be less sensitive. The decrease in sensi- 
tivity with increasing mass can be understood in part by 
noting that for a higher mz^ the cut in Eq. (P^ becomes 
less stringent. For our parameterization the signal rate 
is the same for both Zd masses. Since the cumulative 
background is relatively flat in the region of interest, the 
amount surviving cuts increases. Hence, the significance 
slightly decreases as mz^ increases. Note, however, that 



TABLE II: Luminosities needed for 2a exclusion, 3a observa- 
tion, and 5a discovery after all cuts for both a mz^ = 5 GeV 
and mzj = 10 GeV at y's — 14 TeV with and without K- 
factors. All other parameters are the same as in Eqs. (|12l 
111. 





mza ~ 5 GeV || 




2a (Excl.) 


3a (Obs.) 


5(T (Disc.) 


No ii'-factors 


78 fb^' 


180 fb-" 


490 fb~" 


+A'-factors 


33 fb-' 


75 fb-^ 


210 fb"" 




mz^ = 10 GeV || 




2a (Excl.) 


3a (Obs.) 


5a (Disc.) 


No iC-factors 


100 fb^' 


230 fb"" 


640 fb"" 


-h-ft'-factors 


42fb-i 


95 fb"^ 


260 fb-^ 



lOOOr 



800- 



600- 



B 400 



200 





\ VS = 14 TeV ; 


\ 1 


\ m^ = 5 GeV 




\ m^ = 125 GeV ; 


- 2o\3a 










' '-X^^^.^^ ; 



10' 



10^ 



10'^ 



10'^ 



5^ X Br{Z^-^tr) 



FIG. 6: Luminosities needed for a 2a exclusion (dotted), 3a 
observation (dashed), and 5a discovery as a function of S^ x 
BT{Zd — >■ £^£~) at i/s = 14 TeV. All other parameters are 
the same as in Eqs. p2lll3[) . 



the same S'^ x Br{Zd -^ i~^£~) = 10~^ was used for both 
Zd masses. It is quite plausible, perhaps even likely, that 
S is proportional to mz^- In that case, the 10 GeV signal 
will be enhanced by a factor of 4, requiring 6 — 7 times 
less than the luminosities listed in Table [TTl 

In Fig.[6l the luminosity needed for a 2a exclusion (dot- 
ted), 3ct observation (dashed), and 5a discovery (solid) 
is plotted as a function of S^ x Br {Zd —>■ £'^i~) keep- 
ing mz^ — 5 GeV and all other parameters the same as 
Eqs. ([m [HI). The sensitivity of the LHC dramatically 
decreases as S"^ x Br(Zrf — > £^£^) decreases. 

The effects considered so far have been at leading or- 
der. However, including higher order QCD corrections 
can increase our sensitivity at the LHC. To approximate 
these corrections, some standard iiT-factors ^ are applied. 
The luminosities needed for a 2a exclusion, 3(t observa- 
tion, and 5a discovery after the inclusion of these K- 
factors are shown in Table [Hi As can be seen, the inclu- 
sion of higher order corrections can greatly increase the 
sensitivity of the LHC to these processes, decreasing the 
luminosity needed for an observation by roughly 60%. 

From Table |TT] and Fig. [S] we conclude that with 
~ few X 100 fb~^ the LHC is sensitive to the four lepton 
process H — >■ ZZd for Zd masses in the range 5—10 GeV 
and can probe 6'^ x Br (Zd — )■ £'^£~) down to about 10~^. 
The rescaling of these limits for nonzero kx and kx needs 
to take into account the effects of the cuts. As discussed 
in the next section, the transverse momentum cuts in 
Eq. P^ cut more signal from operators of class (B) than 



'■^ K ^ 2 for gg ^ H a.t NNLO-NNLL in a^ ^]M (N=Next-to; 
LO=Loading Order; LL=Leading Log), K = 1.4 for tf at NNLO- 
NNLL ^, and K = 1.19 for Z at NNLO [H, [Hi. K = 1.18 
for Z7 at NLO in as [MEl, K = 1.62 for ZZ at NLO [13, El, 
K = 0.9 for Zjj at NLO [63 . 
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class (A), due to the transverse polarization of the Z^'s. 
The cuts of Eqs. (|15ll6p and triggers decrease the signal 
rate of operators of class (A) by ~ 50% [see Table [I] and 
operators of class (B) by ^ 65%. These limits can also 
simply be scaled to generalize to the case with different 
branching ratios of Zd into leptons. 



B. Angular Distribution 

As mentioned earlier, interactions mediated by the op- 
erators of class (A) are expected to be dominated by 
the longitudinal polarizations of the Z^,, while the opera- 
tors of class (B) are dominated by its transverse polariza- 
tions. Once a 7? — >■ ZZ^ signal is discovered, determining 
the polarization of the Z^ can allow us to discriminate 
whether the interaction is originating from the operators 
of class (A) or (B). 

A good diagnostic of vector boson polarization is the 
angular distribution of its decay products. In the Zd rest 
frame, we have 



dT{Zd ->■£+£-) 
d cos 9f 



(1 ± cos^ 9e) 



(23) 



where 9e is the lepton {£~) angle with respect to the 
Zd spin-quantization axis, and the upper (lower) sign is 
for a transversely (longitudinally) polarized Zd- Events 
produced via operators of class (A) and class (B) are 
therefore expected to result in different leptonic angular 
distributions of the Zd decay products, allowing for a 
distinction between the two cases. To exploit this effect, 
a spin-quantization axis must be chosen such that the Zd 
is mostly longitudinally or transversely polarized with 
respect to that axis. 

For mzj <C mn — Tnz, the Zd will be produced highly 
boosted and is well approximated as a helicity eigenstate. 
That is, the Zd is mostly in a longitudinally or trans- 
versely polarized state with respect to its moving direc- 
tion. One might expect to be able to use the Zd moving 
direction as the spin-quantization axis and use the angu- 
lar distribution of Eq. (|23p to measure the Zd polariza- 
tion. However, when cosOi ~ ±1, one lepton is moving 
with the direction of motion of the Zd in the lab frame 
and the other lepton in the opposite direction. The boost 
from the Zd rest frame to the lab frame is then against 
the direction of motion of one of the leptons. Hence, this 
configuration results in events with the softest leptons 
and the transverse momentum cut of Eq. p^ eliminates 
the events with cosOe ~ ±1. Since this region is vital in 
distinguishing between our two cases, it will be useful to 
use another angular distribution directly related to that 
ofEq. (f23|. 

In the partonic center of momentum (cm.) frame, the 
Z and Zd momentum and spin must be anti-aligned by 
conservation of momentum and angular momentum, re- 
spectively. As a result, if the Zd is in a helicity eigenstate, 
then in the partonic cm. the Z is also in the same helicity 
eigenstate. The angular distribution of Z decay products 



will be of the same form as Eq. ((23)) . now with the an- 
gle measured in the Z rest frame with respect to the Z 
moving direction in the partonic cm. frame. Since the 
decay products of the Z are typically harder than those 
of the Zd, the pT cuts are not severe and the angular 
distribution is a good diagnostic of the Zd polarization. 

In Fig. [71 we show the simulated angular distribution 
of the lepton identified as originating from the Z de- 
cay for operators of class (A) (solid black) and class (B) 
(solid red/grey), and the expected distributions for lon- 
gitudinally and transversely polarized Zd (dashed black 
and red/grey, respectively). Both operators of class (B), 
Eqs. (|4l5p . result in the same angular distributions as we 
take either k^ 7^ 0, kz = or kz = Q, Hz ^ 0. The CP 
violating effect would show up as an interference effect 
when both operators are present. 

These angular distributions are measured in the re- 
constructed Z rest frame with respect to the Z moving 
direction in the partonic cm. frame. As can be seen 
in Fig. [7] (a), without energy smearing or cuts it is clear 
that, using the above definition, the angular distribution 
of the Z decay products reflect the expected Zd polar- 
ization. That is, operators of class (A) result in events 
that are dominated by longitudinal Z's and Zd's, while 
the operators of class (B) are dominated by transversely 
polarized Z^s and Zd^s. 

Figure [7] (b) illustrates the effects of the cuts on the 
angular distributions. Even after cuts, the two distribu- 
tions are still clearly distinguishable. The small relative 
depletion of events at cos 6'^ ^ ±1 is due to the isolation 
cuts in Eq. (|T6)) . This can be understood by noting that 
for cos 9e ^ ±1 in the Z rest frame, one lepton from the Z 
decay is moving directly against the Z moving direction 
in the partonic cm. frame. Also, since mz ~ rnn, the 
Z is produced nearly at threshold and its decay products 
move back-to-back in the partonic cm. frame. Hence, 
for cos 6*^ ~ ±1, one of the leptons from Z decay is typ- 
ically moving in the Zd direction. In this configuration, 
AR between the Z and Zd decay products is minimized 
and fail the cut in Eq. P^ . 



V. CONCLUSIONS 

Recent data from the LHC appears to have uncovered 
a Higgs scalar associated with the mechanism responsi- 
ble for spontaneous electroweak symmetry breaking. The 
new scalar H may exhibit properties that lead to devia- 
tions from the Standard Model (SM) predictions. Early 
findings suggest inconclusive yet tantalizing hints for such 
deviations in Higgs decays. Regardless of the fate of these 
hints, searching for physics beyond the SM is well moti- 
vated, especially in light of the need to account for the 
cosmic dark matter density in the Universe. It is rea- 
sonable to expect that the "dark" or "hidden" sector of 
particle physics, like its visible counter part, is endowed 
with structure and its own forces. In fact, certain astro- 
physical observations have been interpreted as signals of 




FIG. 7: Solid histograms are event simulation results for angular distributions of the l^ identified as originating from the Z 
in the reconstructed Z rest frame (a) without cuts or smearing and (b) with energy smearing, triggers, and cuts from Eqs. (1151 
- I22|) . The angle is measured with respect to the Z moving direction in the partonic cm. frame. Dashed lines are expected 
distribution for longitudinally (black) and transversely (red/grey) polarized ZdS. 



dark matter that couples to a hidden sector light vector 
boson. One can then ask whether the forces in the dark 
sector could manifest themselves through their interac- 
tions with the Higgs. 

Based on the above motivation, in this work we have 
studied the possibility that a dark vector boson Zd^ in 
the mass range of 5 — 10 GeV, could couple to H via 
mixing or through loop effects. These couplings can then 
be described by two classes of operators, leading to Higgs 
decays into dominantly longitudinal or transverse Zd- An 
interesting typical decay in the first class is -ff — >■ ZZd, 
while the second class of decays could include H -^ XZd, 
with X = Z, Zd, 7. 

We focused on iJ — ^ ZZd, as a representative novel 
decay channel. Using leptonic final states for both Z and 
Zd, we found that the next run of the LHC is capable 
of excluding or detecting such decays for (5^ x Br {Zd — > 
e+i-) - 10-^ and mz^ - 5 - 10 GeV or loop induced 
dimension 5 operators of similar magnitude, with a few 
hundred fb^^ of data. The branching ratios of Zd in 
leptons can be as large as BT{Zd — >■ e+e~) ~ Br(Zd — >■ 
jjL^liT) ~ 0.15. For somewhat lower m^^, larger back- 
grounds will probably require longer running. These 
LHC searches via rare Higgs decays provide a comple- 
mentary approach to low energy experiments that are 
primarily sensitive to dark vector boson masses at or be- 
low the GeV scale. 

Our results also suggest that one could use the angular 
distribution of the leptons from the Z in if — > ZZd to 



probe the underlying micro-physical interaction. Those 
distinct distributions can be used to determine whether 
the Zd emitted in the decay was longitudinal or trans- 
verse, providing a probe of the nature of the Higgs cou- 
pling to the dark sector. 

In the various scenarios we have considered, the Higgs 
boson can also decay into ^Zd and ZdZd- Although we 
have not discussed those modes in detail here, they en- 
tail very distinct signatures [22j that should help sepa- 
rate them from background and ordinary SM Higgs de- 
cays. For example, in the case of very light final state Zd 
bosons, they look like promptly converted e^e~ photons; 
a feature that distinguishes them from ordinary diphoton 
events. Heavier Zd decays could stand out above Dalitz 
and four lepton ordinary decays of the Higgs, depending 
on their abundance in the data. 

We hope that our work encourages further examination 
of the Higgs properties at the LHC, or future facilities, as 
a potential means of peering into the hidden sector and 
shedding light on the nature of dark matter or other as 
yet unknown phenomena. 
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Appendix: Various relevant Higgs decay widths 

1. H ^ ZZa. decay 

The amplitude for H -^ Z{p)Zd{q) is 

^X = tCfizKzz,.^.e*''{p)e*^{q) (A.l) 

where Cf^^ = gmz/ cos 9w is the SM HZZ couphng and 

Kzz^,t,,y = I ezgiiiy H 2" (P • igp^f - Pt^qu) + -^e^.^paP^'q'' ) . (A.2) 

Here, tree-level mixing {ez) is from dimension 3 operators [Eq. (6)] and loop-induced couplings {kz for CP conserving 
part and kz for CP violating part) is from dimension 5 operators [Eqs. (7, 8)]. 
We obtain, after summing over polarizations 

El-^l' = ^C'^zzf (A.3) 

pol 






4 ( S% + 2 1 + 6^^^(P ■<l)+[^] (-l™L + 2(P • <!?) - 2 ( £f ) (rn|„4 - (p • 'Z) 



with p ■ q = h {rn^fi — rn^z ^ "^z ) ^'^^ ^ ^ ZZ^ decay. 
For to| <C D"^ (with D"^ = rn^j — m\), we have 



T{H ^ ZZa) = 4.^^^__ ^ \M\^ (A.4) 



" pol 



(C|zz) £ L_2_2 +3gz^g-T + 47T-T + 4:^-T (A-5) 



with A(a;, y, z) — x'^ + y'^ + z'^ — 2xy — 2yz — 2zx. 

As expected, the longitudinal polarization limit {kztKz — >■ 0) shows an enhancement as mz^ — > 0, in accordance 
with the Goldstone Boson Equivalence Theorem, while the transverse polarization limit {ez — >■ 0) does not. However, 
because ez — Smzj^/mz, the rate is not singular. 

2. H ^ ZdZd decay 

The amplitude for the decay H — > Zd{p)Zd{q) is given by 

tM = iC^ZzKz,z,,,ue*''{p)e*^{q). (A.6) 

When operators of class (A) arise from Z-Zd mass mixing, two insertions of the mixing angle are needed to obtain 
a H-Zd-Zd coupling. This is equivalent to the replacement ez ^ e\. Since operators of class (B) arise from loops, 
the H-Zd-Zd and H-Z-Zd interactions are of the same form in this case. The effective Feynman rule for the H-Zd-Zd 
coupling is of the same form as the right-hand side of Eq. (|A.2|) with ez — > e\ ^'^'^ ^'^ ?fi^ 



pol 



\M\'^{cwizr 



2 / ~ X 2 



.4 (P-I) , o\ I R^2«2d/^ „\ , ( '^Zi\ /_4 , n/^ „\2\ n f ^Za \ /_4 



4|^Td^ + 2)+64;^b-9)+(;^) {m%^+2ip-qy)~2i-^] [m%^~{p-qY 



y I ill y \ IILy I \ IILy 

_(A.7) 
with p ■ q — \ ijn^H ~ 2771^ ) for H -^ ZdZd decay. For the case of scalar mixing in the Higgs sector see Ref. [17| . 
For TO^ <C -D^ (with D^ = rn\j), we have 



r(ff ^ ZdZd) - ^^-^ l^ \M\ (A.8) 

■H pol 

{cTzz) 47-^+34«z.-r + 4.:^-T + 4,i5-T • (A.9) 



?>2TTm% "•""''''' \"^4m|_^ ' "^^'-""ml ' "''"'2m% ' '"^^2^1^ 



11 



3. H —^ ^ Zd decay 



The amplitude for H — > ^i{p)Zd{q) is given by 



(A.IO) 



with the effective Feynman rule of H-^-Zd couphng being the same form as in the right-hand side of Eq. (jA.2p except 
that there is is no ez term here. 
We then have 



pol 



2 u. Z.2\ I ^ 

z 



with p ■ q ^ ^ i^H ~ ^^% ) for _ff — ?► "fZd decay. 
For m| <C -D^ (with D"^ = "^lf)i we obtain 



T[H ^ 7Zd 



47r 



^A(m|^,0,m|J 



TTlr 



(A.ll) 



647r2r 



'ff 



pol 



(C_f/zz) ~3T~r (s + s) 



327r 



mf,m 



(A.12) 
(A.13) 



^'"-2 



As expected, there is no enhancement from the longitudi nal polarization of Zd in this case. 



[1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 

1 (2012) arXiv:1207.7214 [hep-ex]]. 
[2] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 

716, 30 (2012) ;arXiv:1207J23L [hep-ex]]. 
[3] J. R. Ellis, M. K. Gaillard and D. V. Nanopoulos, Nucl. 

Phys. B 106, 292 (1976). 
[4] M. A. Shifman, A. I. Vainshtein, M. B. Voloshin and 

V. I. Zakharov, Sov. J. Nucl. Phys. 30, 711 (1979) [Yad. 

Fiz. 30, 1368 (1979)]. 
[5] "Latest ATLAS studies on Higgs to diboson states," 

talk by Fabrice Hubaut at Rencontres de Moriond 2013, 

EW session, March 2-9; "Study of Higgs Production in 

Bosonic Decay Channels at CMS," talk by Christophe 

Ochando at Rencontres de Moriond 2013, QCD session, 

March 9-16. 
[6] J. Beringer et al. [Particle Data Group Collaboration], 

Phys. Rev. D 86, 010001 (2012). 
[7] C. Boehm and P. Fayet, Nucl. Phys. B 683, 219 (2004) 

[hep-ph/030526i; . 
[8] C. Boehm, D. Hooper, J. Silk, M. Casse and J. Paul, 

Phys. Rev. Lett. 92, 101301 (2004) astro-ph/030968 6 . 
[9] C. Boehm, P. Fayet and J. Silk, Phys. Rev. D 69, 101302 

(2004) ;hep-ph/0311143;. 
[10] P. Jean, J. Knoedlseder, V. Lonjou, M. AUain, J.- 

P. Roques, G. K. Skinner, B. J. Teegarden and G. Ve- 

drenne et al., Astron. Astrophys. 407, L55 (2003) 

[astro-ph/0309484 . 
[11] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer 

and N. Weiner, Phys. Rev. D 79, 015014 (2009) 

,arXiv:0810.071 3 [hep-ph]]. 
[12] J. Chang, J. H. Adams, H. S. Ahn, 



G. L. Bashindzhagyan, M. Christl, O. Ganel, T. G. Guzik 

and J. Isbert et al., Nature 456, 362 (2008). 
[13] O. Adrian! et al. [PAMELA Collaboration], Nature 458, 

607 (2009) arXiv:0810.4995 [astro-ph]]. 
[14] G. W. Bennett et al. [Muon G-2 Collaboration], Phys. 

Rev. D 73, 072003 (2006) hep-ex/0602035 . 
[15] A. Czarnecki and W. J. Marciano, Phys. Rev. D 64, 

013014 (2001) hep-ph/0102122 . 

[16] A. Y. Korchin and V. A. Kovalchuk, arXiv:1303.0365l 

[hep-ph]. 
[17] S. Gopalakrishna, S. Jung and J. D. Wells, Phys. Rev. D 

78, 055002 (2008) |arXiv:0 801.3456 [hep-ph]]. 
[18] H. Davoudiasl, H.-S. Lee and W. J. Marciano, Phys. Rev. 

D 85, 115019 (2012) arXiv: 1203.2947 [hep-ph]]. 
[19] J. M. Cornwall, D. N. Levin and G. Tiktopoulos, Phys. 

Rev. D 10, 1145 (1974) [Erratum-ibid. D 11, 972 (1975)]; 

M. S. Chanowitz and M. K. Gaillard, Nucl. Phys. B 261, 

379 (1985); W. J. Marciano and S. S. D. Willenbrock, 

Phys. Rev. D 37, 2509 (1988). 
[20] H. Davoudiasl, H.-S. Lee and W. J. Marciano, Phys. Rev. 

Lett. 109, 031802 (2012) arXiv:1205.2709 [hep-ph]]. 
[21] K. S. Kum ar, S. Mantry, W. J. Marciano and 

P. A. Souder. larXiv:1 302.6263 , [hep-ex]. 
[22] H. Davoudiasl, H.-S. Lee and W. J. Marciano, Phys. Rev. 

D 86, 095009 (2012) arXiv: 1208.2973 [hep-ph]]. 
[23] M. B. Voloshin, Phys. Rev. D 86, 093016 (2012) 

arXiv: 1208.4303 [hep-ph]]. 
[24] B. Holdom, Phys. Lett. B 166, 196 (1986). 
[25] P. Fayet, Phys. Rev. D 75, 115017 (2007) 

hep-ph/0702176 [HEP-PH]]. 
[26] M. Pospelov, Phys. Rev. D 80, 095002 (2009) 



12 



[27 
[28; 
[29; 

[3o; 
[3i; 

[32; 
[33; 

[34; 

[35] 
[36] 
[37] 
[38] 
[39] 
[40] 
[41] 
[42] 
[43] 

[44] 
[45] 
[46] 

[47] 

[48] 
[49] 

[50] 
[51] 

[52] 



arXiv:0811.1030 [hep-ph]]. 

M. Endo, K. Hamaguchi and G. Mishima, Phys. Rev. D 
86, 095029 (2012) arXiv: 1209.2558 [hep-ph]]. 
J. D. Bjorken, R. Essig, P. Schuster and N. Toro, Phys. 
Rev. D 80, 075018 (2009) arXiv:0906.0580 [hep-ph]]. 
S. Andreas, C. Niebuhr and A. Ringwald, Phys. Rev. D 
86, 095019 (2012) arXiv: 1209.6083 [hep-ph]].' 
D. Babusci et al. [KLOE-2 CoUaboration], 
larXiv: 1210.3927 [hep-ex]. 

R. Meijer Drees et al. [SINDRUM I CoUaboration], Phys. 
Rev. Lett. 68, 3845 (1992). 

S. N. Gninenko, Phys. Rev. D 87, 035030 (2013) 
arXiv:1301.7555 [hep-ph]]. 

P. Adlarson et al. [WASA-at-COSY Collaboration], 
larXiv:1304.067ll [hep-ex] . 

H. Merkel et al. [Al Collaboration], Phys. Rev. Lett. 106, 
251802 (2011) arXiv:1101.4091 [nucl-ex]]. 
S. Abrahamyan et al. [APEX Collaboration], Phys. Rev. 
Lett. 107, 191804 (2011) arXiv;1108.2750 [hep-ex]]. 
R. D. McKeown, AIP Conf. Proc. 1423, 289 (2012) 
[arXiv:110 9.4855 [hep-ex]]. 

S. Dawson and E. Furlan, Phys. Rev. D 86, 015021 (2012) 
larXiv: 1205.4733 [hep-ph]]. 

M. Carena, L Low and C. E. M. Wagner, JHEP 1208, 
060 (2012) arXiv:1206.1082 [hep-ph]]. 
N. Bonne and G. Moreau, Phys. Lett. B 717, 409 (2012) 
|arXiv:1206.3"360l [hep-ph]] . 

H. An, T. Liu and L. -T. Wang, Phys. Rev. D 86, 075030 
(2012) [arXiv:1207.2473l [hep-ph]]. 

A. Joglekar, P. Schwaller and C. E. M. Wagner, JHEP 
1212, 064 (2012) arXiv:1207.4235 [hep-ph]]. 
N. Arkani-Hamed, K. Blum, R. T. D'Agnolo and J. Fan, 
JHEP 1301, 149 (2013) arXi v: 1207.4482. [hep-ph]]. 
L. G. Almeida, E. Bertuzzo, P. A. N. Machado and 
R. Z. Funchal, JHEP 1211, 085 (2012) [arXiv:1207.5254i 
[hep-ph]]. 

J. Kearney, A. Pierce and N. Weiner, Phys. Rev. D 86, 
113005 (2012) arXiv;1207.7062 [hep-ph]]. 
M. A. Ajaib, I. Gogoladze and Q. Shafi, Phys. Rev. D 
86, 095028 (2012) arXiv: 1207.7068, [hep-ph]]. 
S. Dawson, E. Furlan and L Lewis, Phys. Rev. D 87, 

014007 (2013) arXiv:1210.6663" [hep-ph]]. 

H. Davoudiasl, L Lewis and E. Ponton, larXiv:1211.3449l 

[hep-ph]. 

H.-S. Lee and M. Sher. .arXiv: 1303.6653 [hep-ph]. 



K. Agashe, H. Davoudiasl, S. Gopalakrishna, T. Han, G.- 

Y. Huang, G. Pere z, Z.-G . Si and A. Soni, Phys. Rev. D 

76, 115015 (2007) g^v:0709.0007 [hep-ph]]. 

V. Barger, P. Langacker and H.-S. Lee, Phys. Rev. Lett. 

103, 251802 (2009) arXiv:0909.2641 [hep-ph]]. 

J. Alwall, M. Herquet, F. Mahoni, O. Mattelaer and 

T. Stelzer, JHEP 1106, 128 (2011) arXiv: 1106.05221 

[hep-ph]]. 

J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, 



P. M. Nadolsky and W. K. Tung, JHEP 0207, 012 (2002) 
[hep-ph /0201195;. 

[53] N. D. Christensen and C. Duhr, Comput. Phys. Com- 
mun. 180, 1614 (2009) arXiv:0806.4194 [hep-ph]]. 

[54] R. Voss, (ed.) and A. Breskin, (ed.), "The CERN Large 
Hadron Collider, accelerator and experiments," Geneva, 
CERN (2009). 

[55] W.-Y. Keung and W. J. Marciano, Phys. Rev. D 30, 248 
(1984). 

[56] R. V. Harlande r and W. B. Kilg ore, Phys. Rev. Lett. 88, 
201801 (2002) Vp-ph/0201206 . 

[57] C. Anastasiou and K. Melnikov, Nucl. Phys. B 646, 
220 (2002) hcp-ph/0207004 ; V. Ravindran, J. Smith 
and W. L. van Neerven, Nucl. Phys. B 665, 325 (2003) 
[hep-ph/030 2135j ; S. Catani, D. de Florian, M. Grazzini 
and P. Nason, JHEP 0307, 028 (2003) hep-ph/0306211 ; 
V. Ahrens, T. Becher, M. Neubert and L. L. Yang, Phys. 
Rev. D 79, 033013 (2009) arXiv:0808.3008 [hep-ph]]; 
V. Ahrens, T. Becher, M. Neubert and L. L. Yang, Eur. 
Phys. J. C 62, 333 (2009) arXiv:0809.4283 [hep-ph]]. 

[58] M. Czakon, P. Fiedler and A. Mitov, .arXiv: 1303.62541 
[hep-ph]. 

[59] R. Hamberg, W. L. van Neerven and T. Matsuura, 
Nucl. Phys. B 359, 343 (1991) [Erratum-ibid. B 644, 
403 (2002)]; C. Anastasiou, L. J. Dixon, K. Melnikov 
and F. Petriello, Phys. Rev. Lett. 91, 182002 (2003) 
^ep-ph/0306192 ; C. Anastasiou, L. J. Dixon, K. Mel- 
nikov and F. Petriello, Phys. Rev. D 69, 094008 (2004) 
[rep-ph/0312266 ; K. Melnikov and F. Petriello, Phys. 
Rev. D 74, 114017 (2006) hep-ph/0609070 . S. Catani, 
L. Cieri, G. Ferrera, D. de Florian and M. Grazzini, 
Phys. Rev. Lett. 103, 082001 (2009) |arXiv:0903. 21201 
[hep-ph]]. 

[60] J. M. Campbell, R. K. Ellis and C. Williams, JHEP 
1107, 018 (2011) 'arXiv:1105.0Q2Q' [hep-ph]]. 

[61] J. Ohnemus, Phys. Rev. D 47, 940 (1993); U. Baur, 
T. Han and J. Ohnemus, Phys. Rev. D 57, 2823 (1998) 
hep-ph/9710416 ; D. De Florian and A. Signer, Eur. 
Phys. J. C 16, 105 (2000) hep-ph/0002138 . 

[62] J. Ohnemus and J. F. Owens, Phys. Rev. D 43, 3626 
(1991); B. Mele, P. Nason and G. Ridolfi, Nucl. Phys. 
B 357, 409 (1991); T. Matsuura and J. J. van der 
Bij, Z. Phys. C 51, 259 (1991); J. Ohnemus, Phys. 
Rev. D 50, 1931 (1994) hep-ph/9403331 ; L. J. Dixon, 
Z. Kunszt and A. Signer, Nucl. Phys. B 531, 3 (1998) 
^ep-ph/9803250 ; L. J. Dixon, Z. Kunszt and A. Signer, 
Phys. Rev. D 60, 114037 (1999) hep-ph/9907305 . 

[63] J. M. Campbell, R. K. Ellis and D. L. Rainwater, 
Phys. Rev. D 68, 094021 (2003) hep-ph/0308195 ; 

C. F. Berger, Z. Bern, L. J. Dixon, F. Febres Cordero, 

D. Forde, T. Gleisberg, H. Ita and D. A. Kosower et al., 
Phys. Rev. D 82, 074002 (2010) .arXiv: 1004. 1659 [hep- 
ph]]. 



